Exact wave functions for an electron on a graphene triangular quantum dot 
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We generalize the known solution of the Schrodinger equation, describing a particle confined 
to a triangular area, for a triangular graphene quantum dot with armchair-type boundaries. The 
quantization conditions, wave functions, and the eigenenergies are determined analytically. As an 
application, we calculate the corrections to the quantum dot's energy levels due to distortions of the 
carbon-carbon bonds at the edges of the quantum dot. 



I. INTRODUCTION 

Graphene attracts considerable attention due to its 
unusual electronic properties, including: large mean 
free path, "relativistic" dispersion of the low-lying elec- 
tron states, and "valley" degeneracy (see, e.g., reviews 
[1, 2, 3]). These remarkable features suggest that some 
day graphene mesoscopic structures might revolutionize 
nanoscience. Thus, a substantial amount of effort has 
been invested studying graphene nanodevices, such as 
quantum dots (QDs) [4], bilayer structures [5], nanorib- 
bons [6, 7, 8, 9, 10], and other objects: e.g, p — n junc- 
tions, superlattices (both magnetic and non- magnetic), 
and samples with gates [11]. 

In this paper we study graphene QDs, which are in- 
teresting and important nanodevices. A significant char- 
acteristic of a QD is its single-electron spectrum; that 
is, its single-electron wave functions and correspond- 
ing eigenenergies. There is substantial body of litera- 
ture dedicated to investigating the single-electron spec- 
tral properties of graphene QDs using numerical tools 
(see, e.g., Refs. [12, 13, 14]). Instead of using numerical 
approaches, in this paper we obtain an analytical solution 
for a QD shaped as an equilateral triangle (triangular 
QD, or TQD) with armchair-type edges. (Some analyt- 
ical results for a TQD with zigzag edges are reported in 

The basis of our construction is the exact solution of 
the wave equation inside an area shaped like an equi- 
lateral triangle. This solution is considered for differ- 
ent contexts in [16]. The present study is inspired by 
Ref. [13, 14], where the results of Refs. [16] were used 
as a tool of analysis for TQD single-electron numerical 
data. 

Once the wave functions are found, to demonstrate 
their usefulness, we derive corrections to the single- 
electron levels of the TQD due to deformation of the 
carbon-carbon bonds at the edges of the dot. 

This paper is organized as follows. In Sect. II we de- 
scribe the solution of the Schrodinger equation inside a 
triangular well. The necessary basic graphene physics is 



outlined in Sect. III. In Sect. IV analytical expressions 
for the single-electron wave functions for a graphene tri- 
angular quantum dot are found. The properties of these 
wave functions are investigated in Sect. V. The correc- 
tions due to the edge bond deformation are calculated in 
Sect. VI. The obtained results are discussed in Sect. VII. 



II. A QUANTUM PARTICLE INSIDE A 
TRIANGULAR WELL 

Since in this paper we study a triangular graphene 
QD, as a preparatory discussion, let us derive the wave 
function for a quantum particle, confined inside an in- 
finitely deep triangular well. Investigating such a system 
we avoid complications the graphene lattice introduces 
to the problem, yet the most salient features of the wave 
function are brought to light. Thus, we want to solve the 
Schrodinger equation: 

E^/j{x,y) = —ij{x,y), (1) 

with the wave function ^{x,y) vanishing at the bound- 
aries of the equilateral triangle with side L: 

V'(x,0) = 0, (2) 
i}{x,V3x)=0, (3) 
^{x, V3L - VSx) = 0. (4) 

Equation (1) and the boundary conditions Eqs.(2-4) con- 
stitute a well-defined eigenvalue problem. 

As a preliminary step for solving this problem, let us 
ignore Eq. (4) for the time being and construct a wave 
function, which satisfies Eqs. (2) and (3). This amounts 
to solving Eq. (1) inside an infinite sector limited by the 
lines y — and y — \/3x. 

To find the wave function inside the sector we imagine 
that there is an incoming plane wave with wave vector 
ki: 

ipi = exp(-ikir), (5) 

kl = (fca;, ky), (6) 

ky^±V3k^„ ky^O. (7) 




/ 



FIG. 1: Sextet of plane waves which compose the wave func- 
tion ip in Eq. (14). The dashed hnes are ky = ±y/3kx. 



The boundary y = reflects this wave into another plane 
wave 



V'2 = exp(-ik2r), 



with wave vector 



(8) 



(9) 



Now the difference (-01 — ^"2) satisfies Eq. (2). These two 
plane waves are reflected by the boundary y = -s/Sa;, cre- 
ating two additional plane waves V'5,6j whose wave vec- 
tors are: 



k5 = -i(fe 



y j ' 



(10) 
(11) 



These two also experience a reflection at the y = bound- 
ary, inducing two additional plane waves ips^^ with 



k4 — — ~\~ \f^ky^ "s/skx ki 



(12) 
(13) 



Fortunately, when these two undergo reflection at the 
y — ^/3x boundary, no new plane wave appears. The 
sextet of wave vectors (Fig. 1) is closed under reflec- 
tions with respect to the sector's boundaries. A set of 
six plane waves ipa — exp(— ik^r), a = 1 ... 6, is enough 
to describe the wave function inside the sector. The wave 
function in question is: 



(14) 
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FIG. 2: (Color online.) The absolute value \tp{x,y)\ of the 
wave function xp{x,y), Eq. (14), for arbitrary k. The wave 
function vanishes at y = and at y = ±\/^x. 



= E(-l)"^Vo. 



a=l 



By construction it satisfies the boundary conditions 
Eq. (2) and Eq. (3), see Fig. 2. 

Finally, we need to enforce the third boundary condi- 
tion, Eq. (4). On the line 



r = Tq + vs, where 
ro = (L,0), v= i(-l,V3), 



(15) 
(16) 



and s varies from zero to L, our wave function is equal 
to: 

V'(ro + vs) = I [cxp(-ikir) - exp(-ik4r)] (17) 
+ [exp(-zk3r) - exp(-?k2r)] 
-I- [exp(-ik5r) - exp(-ik6r)] \ 

We now group the plane waves exp(— ikir) and 
exp(— ik4r) together because vki = vk4, see Fig. 1 and 
definitions Eqs. (6) and (13). For the same reason we 
cluster exp(— ikar) with exp(— zk2r), and exp(— zksr) is 
grouped with exp(— zkgr). As a result, the value of on 
the boundary can be expressed as: 



0(ro + vs) = Acxp 



-B exp 



^ {kx + VSfcj 



-^ykx "^/iky 



C exp(— ifcjjs) 



(18) 



The coefficients are: 

A = exp(-iA:^L) - exp[i(fca; -I- \/3fcy)L/2], (19) 
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B = exp[i(fca; — V3ky)L/2] — exp{—ikxL), (20) 
C = exp[iikx + V3ky)L/2] (21) 
— e^p[i{kx — VSky)L/2]. 

Equation (17) vanishes, if A, B, and C are all equal to 
zero. This occurs when the following conditions are met: 



kx 



kii 



_ 27r 
27r 

7ii 



[n — m), 
(n + m). 



(22) 
(23) 



Here n and m are integers. Equations (22) and (23) are 
the quantization conditions for the particle momentum 
due to confinement. The wave function ip^'"^{x^y) with 
momentum ki satisfying these equations is the solution 
of the Schrodinger equation with boundary conditions 
Eqs. (2-4) and the eigenvalue: 



E 



8n^n 
9ML'- 



r(n^ + + nm). 



(24) 



The wave function ip^'"^ vanishes identically, if any of the 
equalities 



n = 0, or TO = 0, or n 



-TO, 



(25) 



holds [for example, if n = —to, then fcj, = 0, ki = k2, 
ks = ke, ks = k4, and both the even and odd terms 
of Eq. (14) cancel each other]. Conditions Eq. (25) are 
equivalent to Eq. (7). 

In section IV we show how to adopt for a 

graphene TQD. 



III. BASIC PHYSICS OF A GRAPHENE SHEET 

For completeness, in this section we quickly remind the 
reader the basic single-electron properties of a graphene 
sheet. Our treatment follows Ref. [1]. The notation in- 
troduced in this section will be used in the rest of the 
paper. 

It is common to describe a graphene sample in terms 
of a tight-binding model on the honeycomb lattice. Such 
lattice can be split into two sublattices, denoted by A 
and B. The Hamiltonian of an electron hopping on the 
graphene sheet is given by: 



H 



-*J2 J2 4cr+5, +H.C., 



(26) 



Re^i=i,2,3 



where 'H.c' stands for 'Hermitian conjugate', R runs 
over sublattice A: 

R = di+ aim + SL2n2 <^ R e ^, (27) 

where the primitive vectors of the honeycomb lattice arc: 

ai = ao(3/2,^/3/2), (28) 
a2 = ao(3/2,-V3/2), (29) 



and ni_2 arc integers. The symbol ao denotes the carbon- 
carbon bond length, which is about 1.4 A. The vectors 
Si {i = 1, 2, 3) connect the nearest neighbours. They are: 

<5i = ao(-l,0), (30) 
(52 = ao{l/2,Vs/2), (31) 
S3 = ao(l/2, -\/3/2). (32) 

The corresponding Schrodinger equation can be written 
as: 

= -ti^g+s, - t ^^+5,+a,, (33) 



i=l,2 
i=l,2 



(34) 



where i/'r (V'r+<5i) denotes the wave function value at 
the site R (at the site R-|- (5i) of sublattice A (sublattice 
B). 

The primitive cell of graphene contains two atoms, one 
at R, another at R + ^i. Therefore, it is convenient to 
define the two-component (spinor) wave function: 



(35) 



By construction, the fmiction \I/r is defined on sublattice 
A, Eq. (27). 

The action of on a plane wave 



*R = *k exp(-ikR) 
can be expressed as: 



-tk 



-f* 







tk = t 



l + 2exp(^ „ 2 



3kxao\ ( Vs. 

I r 1 COS I —^kyttQ 



For every k there are two eigenstates: 

1 



*k± = 



exp(i6ik) = 7^, 
Ik 



with eigenvalues: 



ek± = ±|ik| = ±t^-S + F{-k), 
F{\s) = 4 cos (^kxao^ cos ^^^kyaoj 

+2 cos f ^/SkyttQ j 



(36) 

(37) 
(38) 

(39) 
(40) 

(41) 
(42) 



The states with negative (positive) energy are filled 
(empty) at T = 0. 

The allowed values of k lie within the Brillouin zone 
presented on Fig. 3. The reciprocal lattice is character- 
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FIG. 3: The hexagon shown is the Brillouin zone of graphenc. 
The white polygon is where the allowed wave vectors for the 
triangular quantum dot are located (see Fig. 8 and Fig. 9). 

ized by the following lattice vectors: 

di = (47r/3ao,0), (43) 
d2 = (-27r/3ao, 27r/\/3ao). (44) 

The amplitude tk and energy eu± are invariant under 
shifts over di,2- 

The quantity ek± vanishes at the six cor- 
ners of the Brillouin zone: (0,±47r/(3V3ao)) and 
(±27r/(3ao),±27r/(3V3ao)). These are the locations of 
the famous Dirac cones of graphene. 

IV. SCHRODINGER EQUATION SOLUTION 
FOR AN ELECTRON ON A TRIANGULAR 
GRAPHENE DOT 

In this section we find the solution of Eqs. (33) and 
(34) for a graphene TQD. 

The basic object of study here, a TQD with armchair 
edges, is depicted in Fig. 4. The carbon atoms are shown 
as black circles, the covalent bonds are solid lines con- 
necting the atoms. The lateral size of the TQD is Lq. It 
is a multiple of 3ao : 

Lo = SNoao, (45) 

where A^o is an integer. The dot in Fig. 4 is characterized 
by No = 3. 




12a 

FIG. 4: Triangular graphene quantum dot with armchair 
edges (here, Nh = 3 and = 4). The solid linos rcproscut 
covalent bonds between neighboring carbon atoms (black cir- 
cles). The thick lines at the edges represent deformed bonds, 
whose effect on the spectrum is studied in section VI. The 
dashed lines represent fictitious bonds connecting real car- 
bon atoms and auxiliary atoms. The latter arc represented 
by hatched circles. Dotted lines correspond to the effective 
boundaries of the triangular dot. 

The total number of carbon atoms in the dot Na is: 
Na = 3No{No + 1). (46) 

This formula can be derived if one splits the dot into 
No{No + l)/2 aromatic rings, with six atoms each (see 
Fig. 5). 

The atoms at the edges of the dot are special for they 
have only two nearest neighbors, unlike atoms in the 
"bulk" of the dot, which have three neighbors. As a 
result, the Schrodinger equations, Eq. (33) and Eq. (34), 
for the atoms at the edges have to be modified. It is not 
always convenient to work with such formalism. A sim- 
pler approach is used in Ref. [9, 10]. In those works it 
is pointed out that one may add an extra row of carbon 
atoms at the armchair edges ('auxiliary' atoms, shown 
as hatched circles in Fig. 4) and demand the wave func- 
tion to vanish on these 'atoms'. Then for a physical (not 
'auxiliary') atom at the edge we do not have to amend 
Eqs. (33) and (34) explicitly. Indeed, the absent neigh- 
bor (now represented by the 'auxiliary' atom) does not 
contribute to these equations, since the wave function 
vanishes on the 'auxiliary' atoms. 

The addition of the extra row of 'auxiliary' atoms 
slightly increases the effective size of the dot. It is helpful 
to introduce the following notation: 

L = Lo + 3ao = 3Nao, where N = Nq + I. (47) 

Although L and N are trivially related to Lq and Nq, it 

is convenient to define these quantities explicitly for they 
are heavily used in the calculations below. 
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FIG. 5: A triangular graphenc dot can be split into (No + 
l)No/2 aromatic rings, where the integer A^o is proportional 
to Lo, the dot size: Lo = 3aoNo. For the dot shown on the 
figure, No = 3. Thus, the triangular dot consists of six rings. 



Consider now the wave function: 



*R± = X^(-l)"+^*k„±exp(-zk„R), 



(48) 



where are members of a sextet. They are given in 
section II. Observe now that: 



£ki± = ek2± = ek3± = £k4± = £k5± — Skei- 



(49) 



This is a consequence of the graphene lattice symmetry. 
Therefore, the spinor ^'r± is a solution of Eqs. (33) and 
(34) with eigenvalue £ki±- 

Further, the upper component of ^'r± coincides with 
tpCR), Eq. (14). Thus, if ki satisfies Eqs. (22) and (23), 
then tpR± complies with the boundary condition. We re- 
mind the reader that the zero boundary conditions must 
be met at the effective edges of the TQD (on the 'auxil- 
iary' atoms). 

The lower component of ^'r± requires a more tedious 
consideration. It is equal to: 



T 5^(-l)"+i exp(-i0k J exp(-ik„R). (50) 



As Eq. (35) specifies, the argument of tp'^ is not R, which 
belongs to sublattice A, but rather the sum R + (5i , 

which belongs to sublattice B (recall that R is the co- 
ordinate of the two-atom imit cell, while R + is the 
physical location of the atom on sublattice B). We must 
keep this in mind when formulating the following bound- 
ary conditions for ■0'^: 



V'r-f5i± = 0' ifR + (5i = us, 

^n+s,± = 0, if R + 5i = ro + vs', 
where s' = {31 + l)ao. 



(52) 

where u = ^{l,V3),s= {31 - l)ao;(53) 

(54) 
(55) 



Here I is an integer; fq and v are defined in Eq. (16). 

The first condition, Eq. (51), is fulfilled automatically. 
Indeed, it is easy to check that exp(— i6'k) is independent 
of the sign of ky. Thus 

exp(— i^ki) exp(— zkiR) — exp(— z^ks) exp(— ik2R) (56) 
= — 2iexp(— i^ki) exp(— ifca;a;) sm{kyy), 

vanishes, when y = 0. The same holds true for the sum 
of the third and sixth terms, as well as for the sum of the 
fourth and fifth terms. 

Let us now show that Eq. (52) is valid. It is convenient 
to rewrite Eq. (50) as: 

K+5,± = T 5^(-l)"+' exp(-i^k„ - ik^aao) (57) 

a 

X exp [—ika{R + ^i)] • 
When R -h 5i = us, we have: 



V'r-k5i± = Tj^exp 
-l-Bexp 



-^{kx '\/3A^2/)s 



^{kx + V3ky)s 
+ C exp(iA;a:s)|, 



(58) 



where 



A = exp(-i6iki - ikxiao) - exp(-i6'k6 - ikxeao), (59) 
B = exp(-i^k5 - ikxsdo) - exp{-iOu2 - ikx2ao), (60) 
C = exp(-i^k3 - ikxsao) - exp{-idu4 - ikxAao). (61) 

Note that the equation for A involves two wave vectors: 
ki and kg. They enter together because uki = ukg. 
For the same reason the vectors k2 and ks appear in the 
equation for B, and the vectors k3 and k4 are part of the 
equation for C. A similar structure was already observed 
above, see the discussion after Eq. (17). 
It is easy to check that 



exp(-i6'k6) = exp ^-^fc^ao + ^fcj^ao j 

X exp(-i0ki)- 



(62) 



To prove this identity one has to use Eqs. (38) and (40), 
and the fact that |tk| = £k+ is the same for all members 
of the sextet, see Eq. (49). Consequently 



V'R+5,± = 0, if 2/ = 0; 



(51) 



1 — exp 



A = exp(-i6'ki - ikxitto) 
. (3, V3, , , ^ 



(63) 
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Using the definitions of ki and kg we can write: 
_ 3 V3 



(64) 



Thus, the argument of the exponential in Eq. (63) van- 
ishes, and the coefficient A vanishes as a result. In a 
similar fashion, it is possible to prove that B and C are 
equal to zero. 

Lastly, we need to demonstrate that satisfies 
Eq. (54). When R + di = ro + vs', we use Eq. (57) 
to obtain: 



+Bexp 



-ikx + V3ky)s' 



+ C exp{-ikxs')^ , 



(65) 



with the coefficients: 

A = exp[-«6'ki - «fca;(L + flo)] 

- cxp[-«6'k4 + ^(fc^ + V3A;j^)(L-|-ao)], 
B — — exp[—i0]f^2 — ikx{L + Go)] 

+ exp[-i6'k3 + ^{kx - V3ky){L + ao)], 

C = expl-iO),^ + ^{kx + V3ky){L + ao)] 

- exp[-z6'k6 + '^{kx- V3ky){L + ao)]. 



(66) 
(67) 

(68) 



Using the relation: 



3z, 



exp(-z^k4) = exp ^-—kxtto ^kyaoj 

X exp(-i6'ki), 



(69) 



which is similar to Eq. (62) and is derived analogously, 
we show that: 



A = exp[—i9ii^ — ikx{L + ao)] 

a/3, 



(70) 



X 1 1 - exp 
Prom here we obtain: 



iL i ^kx 



A = L (^kx + ^ky 



2-Kn. 



(71) 



For the coefficient B, the following expression holds: 



B = exp[-i^k2 - ikx{L + ao)] 



(72) 



iL ( — fcr 

I 2 



V3, 



X |-l + exp 

Deriving Eq. (72), we use the relation: 

exp(— i^kg) = exp 

X exp(-i6ik2). 



3z , y/li, ^ 

--^kxtto H — ^kyUo 



(73) 



Coefiicient B vanishes when 



L (^kx - -Y'^yj ^ -27rTO. 



(74) 



Coefficient C vanishes automatically, when both Eq. (71) 
and Eq. (74) hold. 

Combining Eq. (71) and Eq. (74) we derive the quan- 
tization condition: 



ki = k"'"*, 
where k"'™ = nKi + mKa, 



Ki,2 



± 



2tt 



2tt 



9Nao ' sVSNao 



(75) 
(76) 

(77) 



which is equivalent to Eqs. (22) and (23) with L given by 
Eq. (47). The symbol ^J^'" is used below to denote the 
wave function with ki = k"'"*. 

Thus, we demonstrate that the wave function Eq. (48) 
with momentum quantized according to Eq. (75) satisfies 
the Schrodinger equations (33) and (34) on a TQD with 
armchair edges. 



V. PROPERTIES OF THE WAVE FUNCTION 

In this section we study the most elementary properties 
of the Schrodinger equation solution ^f^™. 



Eigenenergy 



The eigenenergy corresponding to our solution is equal 



to: 



£n,m± = eki± = ±t <^ 3 -I- 2 COS 



/27rn 
V3iV 



+2 COS 



27rm 



2 COS 



27r(n -I- m) 



3N 



(78) 



1/2 



The eigenenergy remains unchanged, if n and m are 
switched. Therefore, if the wave functions and 
i^m,n linearly independent, the corresponding states 
are degenerate. 

Let us define n and m as: 



n = N + n, 
m = N — fh. 



(79) 
(80) 



If iV > 1, |n| < iV, and \m\ < N, then ^/3N^ « 3N, 
and we may expand Eq. (78) in orders of h/N and rh/N: 



The latter formula is reported in Refs. [13, 14]. 



(81) 
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An interesting phenomenon occurs in TQDs with even 
iV. For such object, consider the quantum states ^'^!^' 
with 



3Ar 



n ■■ 



— m. 



In this case: 



(82) 



(83) 



independent of m. That is, for an even-iV dot the en- 
ergy level at ±t is very degenerate. Although, here we 
do not investigate this feature in detail, it seems to be 
an accidental degeneracy, which is lifted if one includes 
longer-range hopping in the Hamiltonian. 




B. Symmetry of the wave function 

The geometrical symmetry group G of a TQD consists 
of ±27r/3 rotations about the center of the dot and re- 
flections with respect to three bisectors. Such group is 
isomorphic to C^v symmetry group [17]. It has two one- 
dimensional irreducible representations, Ai and A^; and 
one two-dimensional irreducible representation E. 

The representation Ai is trivial: it maps all the group 
elements on 1; ^2 maps all rotations on 1 and all reflec- 
tions on —1. The representation E maps a rotation (re- 
flection) on a 2x2 orthogonal matrix performing a rota- 
tion (reflection) of the two-dimensional Euclidean space. 



1. Rotation 

In order to sec which cigcnfimction corresponds to 
which representation, let us perform a 27r/3 rotation over 
the center of the TQD. 

Technically, it is more convenient to split such trans- 
formation into two consecutive steps: (i) C/27r/3 ~ a ro- 
tation about the origin over the angle 27r/3 (such rota- 
tion does not preserve the location of the dot), followed 
by (ii) a shift over L: x x + L, which restores the 
TQD into its position prior to step (i) (see Fig. 6). After 
step (i) the plane wave ^'k^ exp(— zkiR) becomes the 
plane wave ^'ks exp(— zksR); 5'k2 cxp(— zkiR) becomes 
the plane wave 'fke cxp(— iksR), etc. This means that, 
after the rotation, \1>^™ remains unchanged. 

Next we perform step (ii). As a result of such shift all 
exponentials acquire an extra phase factor. For example, 
consider 

exp(— zkiR) exp(— ikiR — ikxiL) (84) 
= exp{—ikxiL) exp(— ikiR). 

The phase factor is: 



exp(— ifcajiL) = exp(— ifcxi) = exp 



.27r ■ 



(85) 



exp (±^) , if n — m = 3p=F 1, 
1, if n — m = 3p, 



FIG. 6: Sequence of steps to rotate the triangular quantum 
dot about its center. First, the dot is rotated around the 
origin (black arrow). After this transformation, the original 
white TQD becomes the gray TQD. Point A is transformed 
into A', point B is mapped on B'. Afterwards, the dot is 
shifted by L (white arrow) to restore the original position. 



where p is an integer. If we investigate other exponentials 
[exp(ikQ,R), a = 2... 6] we would arrive at the same 
expression for the phase factor. 

Thus, upon rotation around the dot's center, the wave 
function v]/^™ acquires the phase multiplier, Eq. (85). 
We can say that the representation to which the wave 
function belongs is fixed by the value of (n — m). When 
the latter is a multiple of 3, the wave function is trans- 
formed according to Ai or A2. Otherwise, it is part of 
the two-dimensional representation E. 



2. Reflection 

Next, we study how xj/^™ is transformed under reflec- 
tion. Another two-stage process is executed: (a) reflec- 
tion with respect to the x = line followed by (b) shift 
X — > X + L, which restores the original position of the 
TQD. This sequence reflects the dot with respect to its 
vertical bisector. 

After step (a) the wave function becomes 



R 



UR-Si 



(86) 



where ax is the Pauli matrix, and U is the reflection 
transformation matrix: UR = {—x,y). 

Such a complicated transformation law is associated 
with the fact that the reflection x — > —x exchanges 
the sublattices. Thus, the spinor components must be 
switched. This is why we multiply 5* by ax- In addi- 
tion, J7R G B, while the spinor wave function should be 
defined on sublattice A, see Eqs. (27) and (35). Simple 
geometrical considerations show that the unit cell, whose 
location is given by R, is reflected on the cell (t/R — ^i). 
Keeping the above in mind, one can derive Eq. (86). 
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When is a plane wave, Eq. (86) becomes: 

*k± exp(-ikR) (87) 
— > =F exp(— + ik5i)^'[/k± exp(— ikC/R) 
= =F exp(-i0k - ikxaQ)^uu± exp[-i(L''k)R]. 

This equation demonstrates that a plane wave with wave 
vector k is mapped on a plane wave with wave vector 
f/k, multiplied by a phase factor /k = exp{—i9k — ikxao), 
which can be expressed as: 



i 3 i 2 

y exp (ikSa) = > exp 

^+ a=l ^+ s=0 



ik{U2^y5i (88) 



This equation may be proven with the help of Eqs. (38) 
and (40). 

The function / has three important properties: 



if k' 



U2.k 
3 



/k' = 
l/kl 



/k, 
= 1, 



(89) 
(90) 
(91) 



The first property is a simple consequence of Eq. 
while the two others follow from the definition of /. 

Using Eqs. (89) and (91), one demonstrates that, for 
all plane waves in the sextet, the phase factors /k„ are 
identical. It is easy to prove that the transformation law 
for our wave function becomes: 



fn,m — /k^'^ — 



R± 

exp 



-\-Jn,m^Jl± ) 



2m 

"m 



(n — m) 



(92) 
(93) 



1 + expl— nj+exp(-— m 



Equation (92) shows how our wave function is trans- 
formed after step (a) of our two-step process. 

Note that the wave function \1>™'" from the right-hand 
side of Eq. (92) transforms as: 



-rjm.n ^ 



n,m ^ Rib * 



(94) 



Therefore, ^'^!^, subjected to two reflection transforma- 
tions, remains unchanged, as it should be. 

The step (b) is identical to step (ii), see Eq. (84) and 
Eq. (85). Consequently, when the TQD is subjected to 
the reflection about its bisector, the wave function trans- 
forms as follows: 



R± 



~Ffn,m 6xp 



27ri(n — m) 



R± 



(95) 



3. One-dimensional irreducible representations Ai and A2 

At this point we can explicitly construct the wave func- 
tions corresponding to the representations Ai and A2. 



Recall that a wave function belongs to a one-dimensional 
representation {Ai or A2) only when n — m = 3p. As- 
suming this relation, consider the sum: 



^d±a = *R± + c^./' 



(96) 



where a = ±1. Upon reflection, this wave function trans- 
forms as [see Eq. (95)]: 



^R±cr 

Therefore: 



T/n,m*R± T = To-*R±<,- (97) 



Ai : 



R+<7' 



if cr = -1, 



ifa = -l. 



(98) 
(99) 



Since both ^'r!^ and ^^'J^ have identical eigenenergies 
£n,m±j their linear combination ^I/r!^.? ^^^'^ corresponds 

to Sn,m±- 



C. Normalization of the wave function 

In order to calculate matrix elements with the help of 
our wave function, it has to be normalized. Namely, it is 
necessary to flnd the coefficient 7 such that: 



n+Si 



1, (100) 



iRGTQD 



R-F<5ieTQD 



where the summation is performed over the TQD atoms. 

To find 7 we use the following trick. Let us now con- 
sider a large lattice jC, whose linear size is much larger 
than 3A'', the size of our TQD. Consider, further, a spinor 
wave function ^p"'™ on such a lattice (see Fig. 7). This 
wave function vanishes on certain sites of the lattice, 
splitting the whole area into ridot triangular dots. Clearly, 
the sites where the wave function vanishes correspond to 
the auxiliary atoms. 

Using the methods of subsection V B it is possible to 
prove that the wave functions on any two TQD of Fig. 7 
are connected by a unitary transformation. Therefore, 
the summation in Eq. (100), performed over any TQD of 
Fig. 7, gives unity. Thus, the summation over the entire 
lattice C gives us the number of the dots: 



VRG£ 



A\2 



TL+Si I 



= ridot. (101) 



On the other hand, the expression in the round brack- 
ets is equal to 6Nc, where Nc is the number of atoms 
in C. The factor of 6 appears because our wave func- 
tion is composed of six different plane waves. This is the 
advantage of introducing a large lattice: we know that, 
when translational invariance is restored, the interference 
between different plane waves of the sextet is negligible; 
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FIG. 7: (Color online.) The absolute value of the wave func- 
tion vj;"''" on a large graphene lattice. The wave function 
vanishes on the black sites (blue sites when the figure is in 
color), which are the auxiliary atoms. The lines of the auxil- 
iary atoms split the whole lattice in ndot triangular dots. Note 
that a given auxiliary atom is shared by two dots. The bar 
on the right shows the correspondence between the dot color 
and the wave function value, with blue equal to zero and red 
equal to 3.5. 



therefore, each plane wave contributes individually to the 
wave function norm, and no cross-term needs to be cal- 
culated. Thus: 



67 Nc = ndot- 



(102) 



There are Na physical atoms and 3A^ auxiliary atoms per 
TQD on C (there are 6A^ auxiliary atoms surrounding one 
TQD, yet this amount has to be divided by two, since any 
auxiliary atom is shared by two adjacent dots). In total, 
there are 3A^^ lattice sites per TQD. Therefore, we have 



"-dot 



Combining the last two equations we derive: 

1 

7 = 



3V2N' 



(103) 



(104) 



D. Single-electron state labelling 

It appears that for a pair of integer numbers, n and to, 
there is a unique single-electron state. This statement is 
incorrect: not every choice of 71, to is allowed (for exam- 
ple, if n = TO = 0, then the corresponding wave function 
is exactly zero), and not every wave function is unique 
(for example, if we rotate ki by 27r/3 we recover the same 
state). 



It is necessary to introduce a scheme that uniquely 
labels every and any quantum state. The most natural 
way of devising such a scheme is to describe the allowed 
values of k"'™, or, equivalently, of n and to. 

Specifying the allowed k"'™, it quickly becomes obvi- 
ous that the symmetric properties of the sextet are im- 
portant. Therefore, it is convenient to define the sextet's 
symmetry group G. It is isomorphic to C'^y : it consists of 
±27r/3 rotations around the origin and reflections about 
lines ky — 0, ky — zt-v/Sfci- Although, G is isomorphic to 
the TQD's geometrical symmetry group G, they are not 
identical: the reflection axes of G do not coincide with 
those of G. 

Developing this labelling system, one has to abide by 
the following restrictions: (i) if k" ^™ = t/k"'™, where 
L/ e G, then there exists a real number 0, such that 
^n'.m' ^ gi^^n.m. j^^.m ^^^^ within the graphene 

Brillouin zone; (iii) any vector k"'™, such that fc^'™ — 0, 

or fcy^™ = ±'\/3fc"'"', is disallowed: in this case the corre- 
sponding wave function vanishes identically [see discus- 
sion after Eq. (25)]; (iv) there is no state when k"'™ = 
and when k"'™ is the location of the Dirac cone's apex; 
(v) if 



C/k" 



d, [/ e G, 



(105) 



where d is the reciprocal lattice vector, then ^P" ^™ = 

Keeping these conditions in mind let us consider the 
following values for n and to: 



n > 1, 

TO > 1, 

k"'™ € B.Z., 
k"'™ ^ (0,4V(3\/3ao)). 



(106) 
(107) 
(108) 
(109) 



Here 'B.Z.' stands for 'Brillouin zone'. The allowed vec- 
tors k"'™ lie within the white polygon of Fig. 3. For 

= 5 these vectors are shown in Fig. 8. 

Observe that the condition (i) is met: indeed, any two 
allowed vectors, k"'™ and k" '™ , k" ^™ 7^ j^n,™^ cannot 
be connected by G transformations. Conditions (ii)-(iv) 
are explicitly satisfied. 

As for condition (v), it is necessary to realize that it 
is relevant only if k"'™, k" '™ lie on the zone's bound- 
ary. Otherwise, either k"'™, or k" is outside of the 
zone. One can demonstrate that, to satisfy Eq. (105), 
the equality k";^-"" = -/cj?''™' must hold. Thus, in Fig. 8, 
points A and B (open circles) correspond to the identical 
state. The same is true about G and D. 

Finally, we want to count the total number of allowed 
states. It is convenient to group the TQD states as shown 
in Fig. 9. That way they form an arithmetic progression: 
3 states in the first group, 6 states in the second group 
(five states inside the Brillouin zone and one state at the 
zone's boundary), 9 states in the third group, etc. There 
are (A'^ — 1) terms in this progression. The sum of all 
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FIG. 8: The allowed values of k"'™ occupy the sector 
\/3|fea:| < ky of the graphene Brillouin zone (this sector is 
drawn in white in Fig. 3). Every filled circle represents a 
state. Points A and B (open circles near the top) correspond 
to the same state. This is also true for C and D, see subsec- 
tion VD. The thick solid line at the top of the figure is the 
Brillouin zone boundary. 

terms, from the first to the (A'^ — l)th is equal to: 

^Ar(iV-l). (110) 

Since for every k"-" there are two states, and 

the above vahu^ has to be doubled. Therefore, the total 

number of states is equal to 

Nst^tes = 5N{N - 1). (Ill) 

We can see that A^states = -^a- This means that our 
labeUing scheme is exhaustive; that is, there are no states 
unaccounted by it. 

VI. CORRECTIONS DUE TO EDGE BOND 
DEFORMATIONS 

In this section we apply the solution of the Schrodingcr 
equation for a TQD to calculate the correction to the 
single-electron levels due to the deformation of the 
carbon-carbon bonds at the edges of the TQD. 

The edge bonds deformation is known to appear at the 
edges of graphene nanoribbons [7, 8, 9]. The deformation 
is not specific to nanoribbons. Rather, it is a response 
of a carbon-carbon bond to an atypical location (in this 
case, at the edge versus bulk). Thus, it is Ukely that such 
deformation would be present at the edges on a TQD, 
should this device be realized experimentally. 

Our previous calculations completely disregard the 
edge deformation. Fortunately, since the deformation is 




FIG. 9: Counting the quantum states of a TQD. The thick 
solid V-shaped lines show how we group our states to form an 
arithmetic progression. Numbers from 1 to 4 enumerate terms 
of the progression. Note that each pair of open circles (points 
A and B; C and D) counts as one state, see subsection VD. 

weak and since the number of deformed bonds is much 
smaller than the number of undeformed bonds in a suffi- 
ciently large TQD, such modification of the original prob- 
lem can be accounted within the framework of perturba- 
tion theory. Below we show how the deformation of the 
edge bonds affects the single-electron eigenenergies. 

At the Hamiltonian level, we now assume that the 
hopping amplitude across the deformed bond deviates 
from t [8]: 

td = t + St. (112) 

The locations of the deformed bonds are shown in Fig. 4 
by thick solid lines. 

The Hamiltonian due to edge deformations is: 

SH = (ffflowor edge + <5ffloft edge + <5ffri ght cdgci 

(113) 

where the three terms on the right-hand side of the equa- 
tion correspond to the three edges of the dot. 

Let us first discuss the effect due to (fffiowcr edge- The 
deformed bonds at the lower edge connect two atoms 
within the same primitive cell. These cells' positions are 
[see Eq. (27)]: 

R/ ==ao(l/2 + 3Z,V3/2), (114) 
1<1<N-1, integer. (115) 

The matrix element between two arbitrary states ^'r+ 
and $R,+ is equal to: 

N-l 

M = ($+|dffwredge|*+) = $^,+(7,'J'r,4116) 

1=1 
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When $R+ = ^'r+, the matrix element is: 



a,a'=l 



M = St-f^ J2 (-1)"+"' [exp(-i0kj +exp(i0k„,)] (117) 
X ^ exp[-i(k„ - k„')Ri]. 



We can evaluate the sum over I: 

N-l 



^ expH(k„ - k„ORz] (118) 



1=1 



= exp 



N-l 



^ ] exp[ Si{kxa kxa' )aol]. 



1=1 

The sum of the geometric series: 

N-l 



ex.p[-3i{kxa - kxa')(lol] 



1=1 



exp[-3iN{kxa - kxa')ao] - oxp[-3i(fca;a - kxa')a'o\ 



exp[-3i{kxa - kxa')ao] - 1 



(119) 



Depending on a and a', the quantity {kxa — kxa') is equal 
to: 



kxa kx 



_|_ ( 3 L I y^L \ 



(120) 



_ ±VSky. 



Thus, with the help of the condition Eq. (75) we can 
write: 

ro, 

3N{kxa - kxa')ao = s ±TT{n - m)± n{n + m), (121) 
[ ±27r(n + m). 

Therefore, for any a and a' it holds that: 

exp[-3iN{kxa - kxa')ao] = 1 ^ (122) 

N-l 

^ exp[-3i(fca;c« - kxa')aol] (123) 



i=i 



^ r -1, for {kxa - kxa') 0, 
\A''-1, for {kxa - kxa') = 0- 

Expressing the last formula differently, one writes: 



N-l 

^ exp[-3i{kxa - kxa')aol] = -I + NSk^^^k^^,. {124) 
1=1 

The matrix element M can be written as a sum: 

M = Mo + NMi, (125) 



where the Mq term corresponds to —1 in the right-hand 
side of Eq. (124), and NMi terms corresponds to the 
Kronecker delta there. For Mi we obtain: 

6 

Ml = ^ 2 cos6ik, [1 - cos(V3A;j;aao)]- (126) 

a=l 

The first term in the brackets corresponds to summands 
for which a' is such that kxa = kxa' and kya — kya' ■ 
The second term corresponds to summands for which a' 

is StlCll tll&ij kxQ.' cllld kyo^ — — ky(^f . 

To evaluate the sum cos d, it is convenient to use 
Eq. (40) and Eq. (41): 



2 cos 9i,^ = J2 



a=l 



a=l 



(127) 



Since the energy is independent of the index a: e]i^+ = 
ek+, one can write the following expression for this sum: 



At ^ 

' a=l,3,5 



\ / VSky 

a 

1 + 2 COS I — z — I cos 



^2cos0k„ (128) 

ya^O 



a=l 



Substituting the formulas for ki^3^5, Eq. (6), Eq. (12), 
and Eq. (10), into Eq. (128) one obtains: 



2 cos ^ko 



4e 



k+ 



a=l 



t 



(129) 



The calculation of the second term of Eq. (126) is per- 
formed along the same lines. The result is: 



2 cos 0k„ cos(-\/3A:yaao) 

a=l 

t ek+ 



(130) 



cos(3fc2;ao) 



,'3, \ /3V3, \ „ 
+2 cos I -kxao I cos I — ^fcyflo 1 - 3 



Therefore, we can express Mi as follows: 



Ml 



cos(3fcxao) 



(131) 



„ /3, \ /3V3, \ , 
+ 2 cos I -kxao 1 cos I —^kyao I - 3 



where k = {ky,kx), and the function F is defined by 
Eq. (42). 

The evaluation of Mq from Eq. (125) is easy to perform: 
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Mn 



(-l)"+"'exp 



a,a'=l 



')'^0 ^{^ya — kya')ao 



[exp(-i6'k„ ) + exp{i9w^, )] (132) 



= -A7^ ^(-l)"exp --kooaao ^r-kyaao exp(-i6'k„) ^ (-1)" exp -k^oa'ao + -^kya'Uo + C.c. 



c«=l 



a' = l 



r 



where 'C.c' stands for the complex-conjugated terms. 

The sum over a' in Eq. (132) is zero. To prove this let 
us rewrite it: 



6 I (i a/Sz \ 

^ (-1)"' exp -k^a'tto + —kya'tto (133) 



= exp(toouk„/) — exp(iaouka/). 



a'=2,4,6 



a'=l,3,5 



Examining Fig. 1 it becomes obvious that uki = ukg, 
uks = uk4, and uk2 = uk^. This implies that both 
terms on the right-hand side of Eq. (133) are equal, and 
they cancel each other exactly. 

Combining the above results, we write for M: 



M ■ 



2f^NtSt 



[6-F(V3k)]. 



(134) 



This expression gives the matrix clement for the opera- 
tor corresponding to the bond deformations at the lower 
edge of the TQD. The matrix element for the bond de- 
formations at all three edges is equal to 3M: 



3A^£k 

The formula above can be generalized 



[6-F(\/3k)]. (135) 



St 6 - FjVSk) 



to account for the states with negative energies. 

To evaluate first-order corrections to the eigenencrgies 
due to 5H it is necessary to find not only the diagonal el- 
ements ffl!^'™, but the off-diagonal elements, connecting 
the degenerate states, as well. In our case, two wave func- 
tions v]/"''" and vj/™'" correspond to degenerate states, 
unless n — m, OT the vector k"'™ lies on the Brillouin 
zone boundary. However, the element (v[/"''"|5ff |^''"'") 
vanishes. Indeed, SH is invariant under transformations 
from G and, therefore, the matrix element is non-zero 
only if both states transform identically under G. The 
latter condition is never fulfilled, for degenerate wave 
functions either acquire different phase factors upon the 
rotations, Eq. (85), or they transform differently when 
subjected to reflections, Eq. (96) and Eq. (97). 

The above considerations show that the correction to 
the eigenenergies are given by Eq. (136). Let us discuss 
this expression. 



First of all, we notice that [6 - F{^/3k)]/ + F{k) 
is positive. Therefore, the sign of the correction is deter- 
mined by the sign of ±(S. 

Second, since F is an even function of its arguments, 
the degeneracy between 'J'"'"* and remains. 

Third, the larger the dot, the smaller the correction: 

the characteristic energy scale for the correction is St/N, 
which decreases when TV grows. This is natural, since 
the ratio of the deformed bonds (~ N) to the total mun- 
ber of bonds in a TQD (^^ N"^) decreases when the dot 
increases. 



VII. DISCUSSION 



In this paper we find the exact spectrum of a graphene 
TQD with armchair edges. Certain matrix elements are 
evaluated with the help of our wave functions. Thus, 
our solution may be used for perturbation theory calcu- 
lations, e.g., for weak magnetic field, disorder. 

The problem of the electronic properties of graphene 
TQD is addressed numerically in several papers (e.g., [12, 
13, 14]). To show that our analytical approach agrees 
with numerical solutions, we calculated the probability 
density for the states with (a) n = 39, m = 41 (Fig. 10), 
and (b) n = 38, to = 41 (Fig. 11), both for a TQD with 
iV = 41. There are = 4920 atoms in such a TQD. 
The eigenenergies of states (a) and (b) are close to zero. 

These states arc chosen here because their probabil- 
ity distributions are mapped in Ref. [14]. Comparing 
Fig. 4(a) of the latter reference and our Fig. 10 we see 
that the probability density distributions are similar. 
The same is true about Fig. 4(c) of Ref. [14] and our 
Fig. 11. 

To conclude, generalizing the existing solution, we find 
the exact wave functions and eigenenergies for an electron 
inside a graphene TQD. The symmetry properties of our 
wave functions are determined. As an application, the 
corrections to the eigenenergies due to the edge bonds' 
deformations are calculated. We also demonstrate that 
our exact solution is in agreement with previous numer- 
ical work. 
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FIG. 10: (Color online.) Probability density for the state 
with n = 39 and m = 41 for a triangular graphene dot. The 
dot's effective size, L = 3Nao, is fixed by the value of the 
constant A'^ = 41. The total number of atoms in such a dot 
is Na = 4920. The probability density for the same state is 
presented in Fig. 4(a) of Ref. [14]. 
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FIG. 11: (Color online.) Probability density for the state with 
n — 38 and m = 41 for a triangular graphene dot. The dot's 
effective size, L = 3A''ao, is fixed by the value of constant 
= 41. The total number of atoms in such a dot is Na = 
4920. The probability density for the same state is presented 
in Fig. 4(c) of Ref. [14]. 
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